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Resolving the pulsations of subdwarf B stars: 

PG 0154+182, HS 1824+5745, and HS 2151+0857 
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ABSTRACT 

We continue our programme of extended single-site observations of pulsting subd- 
warf B (sdB) stars and present the results of extensive time series photometry to 
resolve the pulsation spectra for use in asteroseismological analyses. PG 0154+182, 
HS 1824+5745, and HS2151+0857 were observed at the MDM Observatory during 
2004 and 2005. Our observations are sufficient to resolve the pulsations of all three 
target stars. We extend the number of known frequencies for PG 0154+182 from one 
to six, confirm that HS 1824+5745 is a mono-periodic pulsator, and extend the num- 
ber of known frequencies to five for HS 2151+0857. We perform standard tests to 
search for multiplet structure, measure amplitude variations as pertains to stochastic 
excitation, and examine the mode density to constrain the mode degree I. 
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1 INTRODUCTION 

Subdwarf B (sdB) stars are thought to be stars with masses 
about 0.5 Mq with thin (<1O _2 M0) hydrogen shells and 
temperatures from 22 000 to 40 000 K (Heber et al. 1984; 
SafFer et al. 1994), making them exceedingly blue. Pulsating 
sdB stars come in two varieties: short period (90 to 600 
seconds) named EC 14026-2647 stars after that prototype, 
officially V361 Hya stars or sdBV stars, with amplitudes 
typically near 1%; and long period (45 minutes to 2 hours) 
named PG 1716+426 stars after that prototype or LPsdBV 
stars, with amplitudes typically <0.1%. For more on these 
stars see Kilkenny (2001) and Green et al. (2003). For this 
work, our interest is the EC 14026 (sdBV) class of pulsators 
as their periods are short, so many pulsation cycles can be 
observed during each run from a single site. 

Pulsating sdB stars potentially allow the opportunity 
to discern their interior structure using asteroseismology, 
obtaining estimates of total mass, luminosity, shell mass, 
radiative levitation, gravitational diffusion, and helium fu- 
sion cross sections. However, to apply the tools of astero- 
seismology, the pulsation frequencies (periods) must first be 
resolved. Variable star discovery surveys seldom resolve or 
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detect the complete set of pulsations. Multisite campaigns, 
because of the complexity of organization, have only ob- 
served a few sdB pulsators. 

Our programme is to resolve poorly-studied sdB 
pulsators, typically from single-site data. This method 
has proven useful for the sdBV stars Feige 48 and 
KPD 2109+4401 (Reed et al. 2004; Zhou et al. 2006). For 
these observations, we obtained most of our data at a modest 
sized telescope (1.3 m) combined with higher signal-to-noise 
(S/N) observations on a 2.4 m telescope. This combination 
allows us to effectively resolve the pulsation spectrum of 
sdBV stars and detect low amplitude pulsations that were 
below the detection limit of the discovery data. Our high 
S/N data from larger telescopes can detect pulsation am- 
plitudes as low as 0.5 milli-modulation amplitudes (mma, 
equivalent to 0.05%) , insuring that we do not miss low ampli- 
tude modes. Our typical temporal resolution is better than 
1.5 /Mz. 

This paper reports the results of our follow-up observa- 
tions on the sdBV stars PG 0154+182 (hereafter PG 0154), 
HS 1824+5745 (hereafter HS 1824), and HS 2151+0857 
(hereafter HS 2151). PG 0154 (£=15.3) was reported as a 
pulsator by Koen et al. (2004; hereafter K04) who detected 
a single frequency at 6090 /iHz. Their observations consisted 
of three < 3 hour observing runs which were obtained un- 
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Table 1. Observations of PG 0154+182 during 2004 on the MDM 
1.3 m telescope 



Run 


UT Start 


Date 


Length 


Integration 




hr:min:sec 


UT 


(hrs) 


(s) 


mdr277 


02:50:30 


06 Oct. 


3.7 


20 


mdr279 


02:50:30 


07 Oct. 


3.8 


20 


mdr281 


02:50:30 


08 Oct. 


3.7 


20 


mdr283 


02:50:30 


09 Oct. 


4.2 


20 


mdr286 


02:50:30 


10 Oct. 


1.2 


20 


mdr288 


02:50:30 


11 Oct. 


4.1 


20 


mdr291 


02:50:30 


12 Oct. 


0.4 


20 


mdr293 


02:50:30 


13 Oct. 


6.6 


20 


mdr296 


02:50:30 


14 Oct. 


0.8 


20 



der less-than ideal conditions and were certainly insufficient 
to resolve the pulsations. From optical spectroscopy and 
2MASS infrared data, they also determined that PG 0154 
has a main sequence G5 - K5 companion. HS 1824 (5=15.6) 
and HS 2151 (B=16.5) were determined to be variable by 
0stensen et al. (2001; hereafter 001) as part of a four star 
discovery paper. (Observations of the other two variables 
will be reported elsewhere.) They detected a single frequency 
in HS 1824 and four in HS 2151. Their observations consisted 
of three « 1 hour runs for HS 1824 and two « 1 hour and 
one 2.5 hour run for HS 2151. As such, the pulsations were 
likely unresolved and ideal for follow-up observations. Both 
HS stars are apparently single with T e g — 33 100 K and 
34 500 K and logg = 6.0 and 6.1 for HS 1824 and HS 2151, 
respectively. 



2 OBSERVATIONS 

Data were obtained at MDM Observatory's 2.4 m and 1.3 m 
telescopes using an Apogee Alta U47+ CCD camera. This 
camera is connected via USB2.0 for high-speed readout: Our 
binned (2 x 2) images had an average dead-time of one sec- 
ond. The observations used a red cut-off filter (BG38) so 
the effective bandpass covers that of B and V filters and 
is essentially that of a blue-sensitive photomultiplier tube. 
Such a setup allows us to maximize light throughput while 
maintaining compatibility with observations obtained with 
photomultiplier photometers. Tables 1 through 3 provide the 
details of our observations including date, start time, inte- 
gration time and run length. 

PG 0154 was observed as a secondary target during our 
fall 2004 KPD 2109 data run (reported in Zhou et al. 2006). 
We obtained short data runs (< 6.6 hr) over nine nights. 
HS 1824 was the subject of a concerted effort, with more 
than 125 hours of data obtained during three different ob- 
serving runs, two at the 1.3 m and the last at the 2.4 m 
telescopes. After the first observing run on HS 1824, we be- 
gan observing HS 2151 as a secondary target. More than 
42 hours of data were obtained of HS 2151 during three 
nights on the 1.3 m and six on the 2.4 m telescopes. Though 
only the observations of HS 1824 conform to our normal 
programme (long time-series observations covering several 
weeks) , we feel that we have successfully resolved all of the 
pulsations. 

Standard procedures of image reduction, including bias 



Table 2. Observations of HS 1824+5745 during 2005 



Run 


UT Start 


Date 


Length 


Int. 


Telescope 




hr:min:scc 


UT 


(hrs) 


(s) 




mdm052505 


04:35:00 


25 May 


6.9 


20 


1.3 m 


mdm052605 


04:12:00 


26 May 


7.4 


20 


1.3 m 


mdm052705 


08:32:00 


27 May 


3.1 


20 


1.3 m 


mdm053005 


04:13:30 


30 May 


7.4 


20 


1.3 m 


mdm053105 


04:43:00 


31 May 


6.9 


20 


1.3 m 


mdm060105 


04:13:30 


01 June 


7.4 


20 


1.3 m 


mdm060305 


03:37:00 


03 June 


7.9 


20 


1.3 m 


mdm060405 


03:48:30 


04 June 


7.7 


20 


1.3 m 


mdm060505 


03:48:30 


05 June 


7.8 


20 


1.3 m 


mdm060605 


03:45:00 


06 June 


7.8 


20 


1.3 m 


mdm061305 


04:17:11 


14 June 


6.6 


20 


1.3 m 


mdm061505 


03:50:00 


15 June 


7.6 


20 


1.3 m 


mdm061605 


03:33:00 


16 June 


8.1 


20 


1.3 m 


hsl8.061705 


03:36:00 


17 June 


3.0 


20 


1.3 m 


hsl8_061905 


03:33:00 


19 June 


1.9 


30 


1.3 m 


hsl8 062005 


03-26 00 


20 lime 


4.4 


25 


1.3 m 


hsl8.062205 


05:40:00 


22 June 


5.9 


25 


1.3 m 


hsl8.070605 


04:17:00 


06 July 


1.2 


7 


2.4 m 


hsl8.070705 


04:14:00 


07 July 


1.5 


7 


2.4 m 


hsl8.070805 


03:21:00 


08 July 


5.4 


7 


2.4 m 


hsl8.070905 


03:18:00 


09 July 


2.2 


7 


2.4 m 


hsl8.071005 


03 15 30 


10 July 


2.1 


7 


2.4 m 


hsl8_071105 


03:20:50 


11 July 


4.1 


7 


2.4 m 


hsl8.071205 


03 23 -30 


12 July 


3.1 




2.4 m 



Table 3. Observations of HS 2151+0857 during 2005 



Run 


UT Start 
hr:min:sec 


Date 
UT 


Length 
(hrs) 


Int. 

(s) 


Telescope 


hs21.061805 


09:10:00 


18 


June 


2.1 


25 


1.3 


m 


hs21.061905 


07:40:00 


19 


June 


3.9 


26 


1.3 


m 


hs21_062005 


08:06:00 


20 


June 


3.5 


28 


1.3 


m 


hs21_070605 


06 01 00 


06 


July 


5.8 


10 


2.4 


m 


hs21_070705 


05:55:50 


07 


July 


5.9 


10 


2.4 


m 


hs21_070805 


08:55:00 


08 


July 


2.9 


10 


2.4 


m 


hs21_070905 


05:43:40 


09 


July 


6.1 


10 


2.4 


m 


hs21_071005 


05:27:00 


10 


July 


6.4 


10 


2.4 


m 


hs21_071205 


06:33:41 


12 


July 


5.2 


10 


2.4 


m 



subtraction, dark current and flat field correction, were fol- 
lowed using IRAF+ packages. Intensities were extracted us- 
ing IRAF aperture photometry with extinction and cloud 
corrections using the normalized intensities of several field 
stars. As sdB stars are substantially hotter than typical field 
stars, differential light curves are not flat due to differential 
atmospheric and colour extinctions. A low-order polynomial 
was fit to remove these trends from the data on a night-by- 
night basis. Finally, the lightcurves are normalized by their 
average flux and centred around zero so the reported dif- 
ferential intensities are Al = (1/(1)) — 1. Amplitudes are 
thus given as milli-modulation amplitudes (mma) with an 



t IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 1. Sample one-hour lightcurves of the three targets. From top to bottom are PG 0154 (obtained on the 1.3 m), HS 1824 (obtained 
on the 2.4 m), and HS 2151 (obtained on the 2.4 m). The dotted line is our least-squares solution of the entire data sets. 



amplitude of 10 mma corresponding to 1.0% or 9.2 milli- 
magnitudes. Sample lightcurves are shown in Figure 1. 



Table 4. Periods, frequencies, and amplitudes for PG 0154. For- 
mal least-squares errors arc in parentheses. 



3 PULSATION ANALYSIS 

PG 0154: The temporal spectrum and window function 
for PG 0154 are plotted in Figure 2. The top panel shows a 
Fourier transform (FT; also known as temporal or pulsation 
spectrum) of the original data, with the two lower panels 
showing the effects of prewhitening by four and six frequen- 
cies, respectively. The inset is the window function, which is 
a single, noise-free sine wave (of arbitrary amplitude) sam- 
pled at the same times as the data. The central peak of the 
window is the input frequency with other peaks indicating 
the aliasing pattern of the data. In order to determine the 
significance of the two low-amplitude peaks, we convolved a 
slightly smoothed window function which included the four 
highest amplitude peaks with the data noise as determined 
using the Breger et al. (1994; hereafter B94) criterion, which 
is 4x the average value of the FT (the dashed line in Fig- 
ure 2). Using such a conservative estimate of the noise, we 
can say with confidence that both low-amplitude peaks are 
real and significantly above the noise. Frequencies, ampli- 
tudes and phases were determined by simultaneously fitting 
a non-linear least-squares solution to the data. Our solution 
for the frequencies and amplitudes is provided in Table 4. No 
other frequencies were detected up to the Nyquist frequency 
near 25000 pHz. 

We readily confirm the frequency detected by K04 and 
detect five additional pulsation frequencies. As the closest 



ID 


Period 


Frequency 


Amplitude 




(s) 


(/iHz) 


(mma) 


/I 


110.9305 (18) 


9014.65 (15) 


1.1 (2) 


/2 


119.5827 (22) 


8362.09 (15) 


1.1 (2) 


/3 


130.0694 (10) 


7688.20 (06) 


2.5 (2) 


/4 


142.2167 (08) 


7031.52 (04) 


3.9 (2) 


/5 


147.3860 (09) 


6784.90 (04) 


3.9 (2) 


/6 


164.2108 (03) 


6089.73 (01) 


9.5 (2) 



two frequencies are separated by nearly 300 /iHz, the pulsa- 
tions in PG 0154 can readily be resolved in about 1 hour of 
observations, a very fortunate circumstance since our data 
runs were short in duration. We also examined individual 
observing runs of sufficient length and readily detected all 
six pulsation frequencies. Though there was some amplitude 
variability (to be discussed in §4.3) which is responsible for 
the peak in the residuals at 6784 /^Hz, the frequencies were 
stable over the course of our observing runs. 

HS 1824: The majority of our 2005 campaign was de- 
voted to confirming that HS 1824 has just a single pulsa- 
tion frequency. The formal least squares solution for the 
frequency, using all of the data is / = 7190.413397 ± 
0.000005 /iHz which is on the low end of 001, but within 
their errors. Using the B94 criterion provides a limit of 
undetected pulsations at 0.48 mma outside of the window 
pattern caused by the single peak. Additional peaks were 
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6000 7000 8000 9000 
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Figure 2. Temporal spectra (FT) for PG 0154 showing the original data (top), and prcwhitcning by 4 and 6 frequencies (middle and 
bottom panels, respectively). Arrows indicate prewhitened frequencies while the dashed line in the middle panel shows the 4<r detection 
limit. Inset is the data window plotted at the same horizontal scale. 



searched for up to the Nyquist frequency near 50000 /iHz, 
but none were detected. However, just because HS 1824 is 
mono-periodic does not make it uninteresting. It has both 
amplitude and phase variations that will be discussed in 
§4. The average amplitude of the pulsation was 2.4 mma 
over the course of our observations, but ranged from 2 to 
6 mma. Such amplitude variation points out a weakness in 
the prewhitening process. Prewhitening algorithms typically 
remove a constant amplitude sinusoidal function from the 
data at the frequency, amplitude, and phase determined via 
non-linear least-squares fitting (or several for multiple pulsa- 
tions). Figure 3 shows how this can affect the prewhitening 
process. A 30 /iHz region centred on the pulsation frequency 
is shown for HS 1824. The top panel shows the original 
FT (the inset shows a 14,000 /iHz span) and the middle 
panel shows the FT prewhitened by one frequency (solid 
line) over-plotted on the window function (dashed line) that 
matches the entire data set. The arrow points to a peak that 
is still above the window function, which could be (mistak- 
enly) deemed intrinsic to the star. The bottom panel shows 
the same prewhitened FT as the middle panel, except that 
the window function is created using the least-squares in- 
formation from individual runs. In this case, there are no 
peaks left above the window and we correctly deduce the 
mono-periodic nature of the pulsations. 

HS 2151: The majority of our data for HS 2151 was 
obtained during our observing run in July, 2005, with just 
three relatively short runs obtained during the late June 
run. As such, we examined these data as three separate sets: 
The June data, the July data, and the combined data. The 
June data alone are not sufficient to resolve the close doublet 
near 7725 /iHz, and the combined data have a slightly more 
complex window function then just the July data. The B94 
noise criterion for the combined and July runs are within 
0.02 mma of each other (0.51 and 0.53 mma, respectively), 
and so the best results were obtained using the July data 
only. The temporal spectrum of HS 2151 is shown in Fig- 
ure 4. The top panel shows the original FT of the July data 
with the next two panels showing the prewhitening sequence 



Table 5. Periods, frequencies, and amplitudes for HS 2151. For- 
mal least-squares errors are in parentheses. 



ID 


Period 


Frequency 


Amplitude 




(■) 


(/iHz) 


(mma) 


/I 


129.4133 (07) 


7727.18 (4) 


3.01 (13) 


/2 


129.4671 (13) 


7723.97 (8) 


1.51 (13) 


/3 


134.6929 (17) 


7424.30 (9) 


1.20 (12) 


/4 


145.7850 (13) 


6859.41 (6) 


1.74 (12) 


/5 


151.1561 (07) 


6615.68 (3) 


3.84 (12) 



of two and five frequencies, respectively. The insets show 
the window function and an enlarged section of the origi- 
nal FT to show the resolved doublet. Arrows indicate the 
prewhitened frequencies and our least-squares solution is 
provided in Table 5. No further frequencies were detected 
up to the Nyquist frequency near 39000 /iHz. 



4 DISCUSSION 

4.1 Comparison with the discovery data 

The goal of these observations was to obtain better data 
than in the discovery papers and to resolve the pulsation 
frequencies for asteroseismic analysis. For the sake of com- 
parison with the discovery data, we calculate the temporal 
resolution as 1/AT with AT being the temporal extent of 
the observations. For the discovery data, we can determine 
temporal resolution from information provided in their pa- 
pers, and estimate "by eye" the discovery data detection 
limit as twice the top of the noise level in their FTs outside 
of the pulsations and windows. 

For PG 0154, we were able to confirm the frequency 
from K04 as well as detect five new frequencies. Our obser- 
vations have a temporal resolution of 1.4 /iHz which is more 
than 4x better than the discovery data and a detection limit 
of 0.76 mma which is more than 6x better than in K04. 
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Figure 3. Temporal spectra (FT) for HS 1824 showing the pitfalls of prcwhitcning. The top panel shows the original FT (the inset 
spans a larger region) and the middle panel shows the original FT prcwhitcned by one frequency (solid line) using a window function fit 
from the entire data set (dashed line). The bottom panel shows the same FT as the middle panel, but the window function incorporates 
the amplitude and phase variations determined from individual nights of data. 




Figure 4. Temporal spectra (FT) for HS 2151 showing the original data (top), and prewhitening by 2 (middle) and 5 frequencies 
(bottom). Arrows indicate prewhitened frequencies. Inset shows the data window and an expanded view of the close pair of frequencies. 



We certainly obtained more and better data than 001 
for HS 1824, but it ended up only confirming the singe fre- 
quency they detected. Our temporal resolution is 0.25 /iHz, 
which is 23 x better than theirs and our detection limit is 
0.48 mma which is w 4x better than 001. 

For HS 2151 we confirm the four previously detected 
frequencies of 001 and detect an additional, closely spaced 
frequency. This detection was possible because of our 11 x 
better temporal resolution of 0.5 /jMz and our detection limit 
was also « 3x better than 001 at 0.53 ^iHz. 

Compared to the discovery data of these three stars, we 
certainly obtained our goal of vastly superior temporal reso- 
lution and improved detection limits. However, it is interest- 
ing to note that both PG 0154 and HS 1824 showed a single 
peak in very limited discovery data, and while PG 0154 paid 
dividends with only relatively few observations, no new pul- 
sations were detected in HS 1824 after a concerted effort. 



4.2 Multiplet constraints on pulsation modes 

For our goal of applying asteroseismological tools to these 
stars, one of the most helpful features would be observational 
constraints on the pulsation modes themselves, rather than 
just recording the pulsation frequencies. Pulsation modes 
can be described by three quantum numbers, n (or k), £, 
and m, and mathematically described by spherical harmon- 
ics. As rotation can break the m degeneracy by separating 
each degree £ into a multiplet of 21+1 components, multiplet 
structure is a very useful tool for observationally constrain- 
ing pulsation degree (see Winget et al. 1991). Such multi- 
plets should be nearly equally spaced in frequency but such 
structure is seldom observed in sdBV stars. However, when 
multiplets are observed, they not only constrain the pulsa- 
tion degree (£) but also the rotation period and inclination 
(Reed et al. 2004; O'Toole, Heber, & Benjamin 2004). 

For PG 0154 there is a common value that emerges with 
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Table 6. Frequency spacings between pulsation frequencies for 
PG 0154. 

Frequency pair Difference 



(MHz) 


/I 


-/2 


652.6 


/2 


-/3 


673.9 


/3 


-/4 


656.7 


/5 


-/6 


695.4 



an average frequency spacing of ss 670/iHz that includes all 
of the observed frequencies. The spacings are provided in 
Table 6 and range from 652 to 695 /iHz. This is less than 
a 5% variation and roughly of the order expected from the- 
ory. If these are rotationally split multiplets, /l through /4 
would necessarily belong to a single degree of I > 2, miss- 
ing at least one component (just one for £ — 2) and /5 and 
/6 would be two parts of an £ > 1 multiplet. As such, all 
of the frequencies in PG 0154 can be explained using one 
£ = 1 and one £ = 2 mode; each missing one of their com- 
ponents. However, rotational splitting of 670 /iHz means a 
rotation period of 1492 s, or about 25 minutes! At a canon- 
ical radius of O.15R,0, this would imply a rotation velocity 
of w 400 km/s; an easily measurable quantity. Spectroscopic 
constraints already exist and from Figure 7 of K04, it is ob- 
vious that PG 0154 is a normal, slowly rotating sdB star. 
Though it could be possible to conceal such a high veloc- 
ity with a low inclination, another argument against this 
scenario is the orbital separation. Using the canonical sdB 
mass of 0.5 Mq and a range of masses appropriate for the 
main sequence companion, the orbital distance between the 
two mass centres is <0.95 R©. The two stars could not exist 
within such close proximity to each other without significant 
effects; at a minimum, the sdB star would have significant 
ellipsiodal variations like those observed for the close binary 
KPD 1930+2752 (Billeres et al. 2000) and for most com- 
panion parameters, the two stars would not fit within the 
orbital separation. 

Since the spacing is so large, it is conceivable that they 
could be successive overtones of the radial index, n. How- 
ever, because /4 and /5 do not share this spacing, they 
would necessarily belong to differing degrees, £, and thoery 
predicts that differing degrees have different overtone spac- 
ing. However, a glance at Figure 4 of Charpinet et al. (2002) 
reminds us that theory expects these to be low overtone 
pulsations, which do not obey the asymptotic relationship 
of even frequency spacing. The figure indicates that succes- 
sive overtones should differ in spacing by at least 500 /iHz, 
ruling out this possibility. 

Another possibility is the "Kawaler scheme" which has 
recently been introduced (Kawaler et al. 2006; Vuckovic et 
al. in preparation). Though it is not compatible with low 
overtone (n) pulsation theory invoked for sdBV stars, it has 
been noticed that for several stars, an improved frequency 
fit can be obtained using an asymptotic-like formula; 

f(hj) = fo + ixS + j x A 

where i has integer values, j is limited to values of —1, 0, and 
1, and up till now S has been a small spacing and A has been 
a large spacing. However, for PG 0154, the large spacing 
occurs too frequently to be restricted to just three values, 
( — 1, 0, +1). For PG 0154, there are not enough frequencies 



Table 7. A "Kawaler scheme" model fit to the observed frequen- 
cies using f = 6349.9 (6103.2), <5 = 669.6 and A = 247.3/iHz. 



Star 


i 


j 


Model 


Difference 


/I 


4 


0(1) 


9028.1 


-13.5 


/2 


3 


0(1) 


8358.6 


3.5 


/3 


2 


(1) 


7689.0 


-0.8 


/4 


1 


(1) 


7019.4 


12.1 


/5 


1 


-1 (0) 


6772.5 


12.4 


/6 





-1 (0) 


6103.4 


-13.1 



to make a unique fit and so Table 7 gives the solutions for two 
possible values of f a and Figure 5 shows the Echelle diagram 
of the fits using a different symbol for each fit. (Fractional 
multiples of the S spacings, i.e. 8/2, are also equally good 
fits but the largest spacing comes directly from the data.) 
For both models, S = 669.6/iHz and A = 247.3/iHz using 
fo = 6349.9/iHz for the squares and / = 6103.2/iHz for 
the triangles. Part of the Kawaler scheme is to choose "the 
(5, A) pair for which we have at least two modes with the 
same value of i but different j" and we have satisfied this 
condition, but without coming to a unique conclusion. Still, 
we can say that the Kawaler scheme fits the data with an 
r.m.s. error of only 0.16% and remains an option while a 
rotational explanation for the frequency splitting is clearly 
ruled out. A Kawaler scheme fit was also obtained switching 
the values for 5 and A, but the differences were substantially 
larger. 

For HS 2151, where we have five frequencies to work 
with, no two frequency spacings are similar. As such, we are 
left with the following possibilities for the five pulsation fre- 
quencies detected: 1) Rotation is sufficiently slow that all m 
values are degenerate; 2) at most one pair share the same n 
and £ values with differing m, which seems very unlikely; 3) 
our line of sight is along the pulsation axis, with sini ~ 0, 
making only m = modes observable because of geometric 
cancellation (Pesnell 1985; Reed, Brondel, & Kawaler 2005); 
or 4) internal rotation is such that rotationally-induced mul- 
tiplets are widely spaced and uneven (Kawaler & Hostler 
2005). 

4.3 Constraints on pulsation degree via mode 
density 

An open question involving sdBV stars is the mode degree 
£ of the pulsations. In resolved sdBV stars, we sometimes 
observe many more pulsation modes than £=0, 1, and 2 can 
provide. Higher £ modes may be needed, but if so they must 
have a larger intrinsic amplitude because of the large de- 
gree of geometric cancellation (Charpinet et al. 2005; Reed, 
Brondel, & Kawaler 2005) . A general guideline would be one 
n order per £ degree per 1000 /iHz. As such, the temporal 
spectrum can accommodate three frequencies per 1000 /iHz 
without the necessity of invoking high-<? values. 

For PG 0154, even ignoring the frequency spacing that 
inter- relates several frequencies, it is still easy to accommo- 
date all the frequencies without the need for £ > 3 modes. 
The pulsation spectrum for HS 2151 is slightly more dense 
and we do not detect any related frequencies. However, it is 
still possible to fit the frequencies within the range observed 
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Figure 5. Echelle diagram for PG 0154 using "Kawaler scheme" model fits. Squares correspond to the model with f 
triangles for the model with f a = 6103.2 fiHz. 



- 6349.9 ^Hz and 



without invoking £ > 3 modes if the 6616 and 6859 /iHz 
frequencies are related to those near 7725 /iHz. As such, we 
can deduce that the pulsations in all three pulsators in this 
paper can be accommodated using I = 0, 1, and 2 modes. 
However, this by no means implies that they are. 



4.4 Amplitude variability 

Another feature we can examine is the amount of amplitude 
variability in the pulsations. If pulsating sdB stars are ob- 
served over an extended time period, it is common to detect 
amplitude variability in many, if not all, of the pulsation 
frequencies. Such variability can occasionally be ascribed to 
beating between pulsations too closely spaced to be resolved 
in any subsets of the data, but often appear in clearly re- 
solved pulsation spectra where mode beating cannot be the 
cause. 

Figure 6 shows the amplitudes of the pulsations of 
PG 0154 for individual nights. The scale is selected (for 
all of the amplitude and phase plots) such that a 10% er- 
rorbar occupies ^10% of the panel. All of the amplitudes 
are essentially constant to within the errors except for /5, 
which triples in amplitude over the course of the run. We 
initially suspected it was caused by an unresolved doublet 
of frequencies, which could still be the case. However, several 
other frequencies show the same general trend (though on 
smaller scales) and the amplitude variations do not appear 
sinusoidal, which would be expected with mode beating. As 
such, we doubt our initial interpretation and rather think 
these variations are intrinsic. However, the only way to be 
sure is to obtain more observations. Still, the net result is 
that five of six frequencies are stable both in amplitude and 
phase (not shown as all of the phases are constant to within 
the errors). 

Of significantly more interest are the amplitudes and 
phases of the single frequency present in HS 1824. These are 
shown in Figure 7 and indicate that they are quite vari- 
able. Between observing days 41 and 45 (JD=2453531.5 
and 2453536.5), the amplitude rises from 0.9 to 4.2 mma 
while the phases change by w 40% between days 42 and 



46. Clearly, HS 1824 is both amplitude and phase variable; 
possibly an indicator of stochastic pulsations (discussed in 
the next section). Note that phases are chosen such that the 
phase is the time of first maximum since JD=2453279.5 for 
PG 0154, JD=2453491.5 for HS 1824 and HS 2151 divided 
by the period. 

Figure 8 shows the amplitudes and phases of the three 
well separated frequencies of HS 2151. The doublet near 
7725 /iHz is not resolvable on a nightly basis. In this case, 
two of the three frequencies (/4 and /5) have amplitude 
variations marginally larger than the errors which appear to 
follow the same trend. One possible explanation would be 
if they are two parts of a multiplet sharing an increase in 
driving power, though such a hypothesis cannot be tested 
using these data. For all three frequencies, the phases vary 
over the duration of our observing campaign. The highest 
amplitude frequency (/5) is generally the most stable; vary- 
ing smoothly by 15% over the first three nights and then 
remaining essentially constant. The other two frequencies 
(/3 and /4) both show an offset between the June and July 
data runs with /4 essentially constant during the individual 
month's data. This points toward a problem with timing 
rather than intrinsic changes of the star. We investigated 
this possibility but ruled it out for the following reasons: 
The highest amplitude frequency does not show these, the 
changes are not the same in period (i.e., the same amount of 
difference in timing), our observations of HS 1824, which oc- 
curred during the same nights do not show this shift (though 
admittedly HS 1824 shows phase variations), and the timing 
setup is the same as we have used during all of our other 
runs^ and is checked multiple times during the night. As 
such, we conclude that these shifts that occurred during the 
15 night interval between observations are intrinsic to the 
star. 



t The acquisition computer's clock is synchronized with time 
servers through Network Time Protocol. 
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Figure 6. Pulsation amplitudes for frequencies detected in PG 0154. Error bars are the formal least-squares errors. 
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Figure 7. Pulsation amplitudes and phases for frequencies detected in HS 1824. Error bars are the formal least-squares errors. 
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Figure 8. Pulsation amplitudes and phases for frequencies detected in HS 2151. Error bars are the formal least-squares errors. 
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Table 8. Mean amplitude, standard deviation, and their ratio for 
readily resolvable frequencies. 



J_ 1 LlCHt, V 

(MHz) 


I A\ 
(mma) 


rr( A\ 
(mma) 


rr( A \/l A\ 




PG 


0154 




9015 


0.97 


0.38 


0.39 ± 0.39 


8362 


1.37 


0.34 


0.25 ± 0.22 


7688 


2.72 


0.36 


0.13 ± 0.12 


7032 


4.07 


0.49 


0.12 ± 0.05 


6785 


3.90 


1.27 


0.33 ± 0.06 


6090 


9.53 


0.48 


0.05 ± 0.03 




HS 


1824 




7190 


3.19 


1.01 


0.31 ± 0.11 




HS 


2151 




7424 


1.24 


0.16 


0.17 ±0.12 


6859 


1.77 


0.40 


0.22 ± 0.10 


6616 


4.18 


0.55 


0.13 ± 0.04 



4.5 The nature of the excitation mechanism 

We can also use the criterion outlined in Christensen- 
Dalsgaard et al (2004; hereafter JCD04) and first applied 
to sdBV stars in Pereira & Lopes (2005) to examine the 
nature of the excitation mechanism. If the frequencies are 
stochastically excited, we can expect to find the standard 
deviation of the amplitudes, o(A), divided by the average 
amplitude (A) to have a value near 0.5 (Eqn. 7 of Pereira & 
Lopes 2005). 

Both parameters and their ratios are given in Table 8 
for all resolved frequencies in our target stars. For PG 0154, 
three of the six frequencies have relatively stable pulsation 
amplitudes, and their ratios reflect this as non-stochastically 
excited pulsations. The two low amplitude frequencies, /l 
and /2, have errors too large to constrain the ratio, which 
for /5 is in between. In this case /5 suffers the same problem 
as it did for mode beating: The time-scale of amplitude vari- 
ation appears longer than our observations. As such, more 
data would be required to discern which interpretation is 
correct for /5: Stochastically driven pulsations, mode beat- 
ing, or variable amplitudes caused by something else. Sur- 
prisingly the ratio for HS 1824 is nearly the same as for 
/5 of PG 0154. However we can suggest that these pulsa- 
tions do not match the JCD04 criteria for stochastic oscil- 
lations as our observations are surely sufficient in duration 
to include several damping lifetimes, yet a(A)/(A) is only 
0.31 ±0.11. As such, the rather large variation in amplitude 
and phase remains unexplained. For the three well-separated 
frequencies of HS 2151, none of them have ratios indicating 
stochastically excited pulsations. These results indicate that 
at most one of our ten resolved frequencies have indications 
of being stochastically excited. This can be compared to 
our previous observations of KPD 2109±4401 where 50% of 
the resolved pulsations indicated stochastic excitation and 
a study of PG 1605±072 where no indications of stochastic 
excitation were found from 11 resolvable frequencies (Pereira 
& Lopes 2005). It is a bit early to interpret the meaning of 
such results and it is not clear what impact the discovery 
of many stochastically driven frequencies (should that oc- 
cur) would mean for the proposed iron driving mechanism 
of Charpinet et al. (2001). 



5 CONCLUSIONS 

From extensive follow-up data acquired at MDM observa- 
tory, we are confident that we have resolved the pulsation 
spectra of three additional pulsating sdB stars. We are able 
to confirm all of the previously observed frequencies as well 
as detect a total of six new frequencies in two of our three 
targets. 

For PG 0154, we have discovered five new frequencies 
with the closest spacing between frequencies at ~ 250 /iHz. 
As such, they are readily resolvable in 1.5 hours of data, 
however the low amplitude of two frequencies requires high 
S/N or longer time-base observations to reduce the noise. We 
notice an enticing common frequency spacing of ss 670 fiRz, 
which could be interpreted as all six frequencies being mem- 
bers of one I = 1 and one I — 2 pulsation multiplet. How- 
ever, this would require a rotation period near 25 minutes 
and a rotation velocity of ~ 440 km/s which is clearly ruled 
out by K04. We applied the Kawaler scheme to the spacings 
using a purely numerical fit with good results. However, as 
PG 0154 only has six frequencies, there is no unique solution 
and so we are not convinced that it should be applied in this 
case. As such, the cause for the large, nearly-equal frequency 
spacing remains a mystery. Three of the six frequencies are 
amplitude-stable, while one shows variations that could be 
indicative of a stochastically excited mode. Of course this 
would not be possible if it is actually an m component of 
the same degree as other, amplitude-stable frequencies. As 
such, it seems we still have some to learn about PG 0154 and 
additional follow-up observations could still be instructive. 

In HS 1824 we confirmed the frequency detected by 
001, but did not uncover any additional frequencies. How- 
ever, we were able to measure significant amplitude and 
phase variations over a long span of time. As such, once 
asteroseismology matures to where it can interpret such fluc- 
tuations, HS 1824 will be an excellent target. But until then, 
all we can readily say is that the variability in the amplitude 
and phase of pulsation seems unrelated to stochastically ex- 
cited oscillations. 

For HS 2151, we confirm the four previously known fre- 
quencies and detect an additional oscillation close to a pre- 
viously observed one. However, none of the frequencies ap- 
pear to have common spacing indicative of multiplet struc- 
ture and their amplitudes appear stable, indicating non- 
stochastically excited pulsations. 

These three pulsating sdB stars show the variety, and 
perhaps the complexity, of oscillations observed in sdB stars. 
Two of the three have several frequencies, while HS 1824 
only had one. Of the ten total frequencies, several show am- 
plitude variations, but only one has sufficient variability to 
even consider stochastic excitation. Several of the phases 
show small variations with those in HS 1824 being much 
larger and only one star shows a common frequency spacing 
that could be used to observationally constrain the pulsation 
degree, £, except it would require a drastically fast rotation 
velocity which has been ruled out by previous spectroscopic 
observations. 
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